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ABSTRACT 

Attempts to consider evolution across space-time singularities often lead to quantum 
systems with time-dependent Hamiltonians developing an isolated singularity as a function 
of time. Examples include matrix theory in certain singular time-dependent backgounds 
and free quantum fields on the two-dimensional compactified Milne universe. Due to the 
presence of the singularities in the time dependence, the conventional quantum-mechanical 
evolution is not well-defined for such systems. We propose a natural way, mathematically 
analogous to renormalization in conventional quantum field theory, to construct unitary 
quantum evolution across the singularity. We carry out this procedure explicitly for free 
fields on the compactified Milne universe and compare our results with the matching con- 
ditions considered in earlier work (which were based on the covering Minkowski space). 



1 Introduction 

Dynamical evolution across space-time singularities is one of the most tantalizing, even if 
speculative, questions in modern theoretical physics. Should our theories point towards a 
beginning of time, it is very natural to ask what came before, and, indeed, whether there 
could be anything before. 

In certain model contexts, quantum evolution across space-time singularities appears 
to be described by time-dependent Hamiltonians developing an isolated singularity as a 
function of time at the moment the system reaches a space-time singularity. It is then 
worthwhile to study such quantum Hamiltonians and establish some general prescriptions 
for using them to constuct a unitary quantum evolution. Needless to say, additional spec- 
ifications are needed in a Schrodinger equation involving this kind of Hamiltonians, on 
account of the singular time dependence. 

One of the simplest examples of such singular time-dependent Hamiltonians in systems 
with space-time singularities is given by a free scalar field on the Milne orbifold (see [TJ [21 
|3l HI [5] and references therein for some recent occurrences of the Milne orbifold in models of 
cosmo logical singularities). We shall give a detailed consideration of this case in section 3. 
Here, it should suffice to say that the square root determinant of the metric of the Milne 
orbifold vanishes as \t\ when t goes to 0. Because of that, the kinetic term in the Lagrangian 
for a free field on the Milne orbifold will have the form \t\(d t 4>) 2 , and the corresponding 
term in the Hamiltonian expressed through the canonical momentum conjugate to 
will have the form 7r|/|t|, which manifestly displays an l/\t\ singularity. The position of 
this singularity in the time dependence coincides with the metric singularity of the Milne 
orbifold. 

While it is well-known that free fields on the Milne orbifold are not a good approximation 
to interacting systems, especially in gravitational theories [6J [7J, analogous singular time 
dependences have recently appeared in other models, which have been the main motivation 
for the present work. For example, 11-dimensional quantum gravity with one compact 
dimension in a certain singular time- dependent background with a light-like isometry is 
conjectured to be described by a time-dependent modification of matrix string theory [SJ, E] ■ 
This model can be recast in the form of a (l+l)-dimensional super- Yang-Mills theory on 
the Milne orbifold. It will thus contain in its Hamiltonian the l/\t\ time dependence typical 
of the general Milne orbifold kinematics. The question of transition through the singularity 
will then amount to defining a quantum system with such singular Hamiltonian. Likewise, 
for the time-dependent matrix models of [10J, which are conjectured to describe quantum 
gravity in non-compact eleven- dimensional time- dependent background with a light-like 
singularity, one obtains a quantum Hamiltonian with a singular time dependence. 

In view of these examples, our present paper will address the question of how one 
should define unitary quantum evolution in the presence of isolated singularities in the 
time dependence of quantum Hamiltonians. Upon giving a general prescription for treating 
such singularities and discussing the ambiguities it incurs, we shall proceed with analyzing 
the simple yet instructive case of a free scalar field on the Milne orbifold. We shall further 
discuss the relation between our prescription and the recipes for quantum evolution of this 
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system previously proposed in the literature (and based on considerations in the covering 
Minkowski space) pH H21 H31 Ej . 



2 Isolated singularities in time-dependent quantum Hamiltoni- 
ans 

Following the general remarks in the introduction, we shall consider a quantum system 
described by the following time-dependent Hamiltonian: 

H(t) = f(t,e)h + H reg (t), (1) 

where H reg (t) is non-singular around t — 0, whereas the numerical function f(t,e) develops 
an isolated singularity at t = when e goes to (e serves as a singularity regularization 
parameter), and h is a time-independent operator. We shall be interested in the evolution 
operator from small negative to small positive time. In this region, we shall assume that 
we can neglect the regular part of the Hamiltonian H reg {t) compared to the singular partQ 
The Schrodinger equation takes the form 

ij t \*) = f(t,e)h\*). (2) 
The solution for the corresponding evolution operator is obviously given by 



U(t, t') = exp 



dtf(t,e)h 



(3) 



When the regularization parameter e is sent to 0, f(t,e) becomes singular and U(t,t') is in 
general not well-defined. 

The goal is then to modify the Hamiltonian locally at t = in such a way that the 
evolution away from t — remains as it was before, but there is a unitary transition 
through t — 0. Of course, a large amount of ambiguity is associated with such a program, 
and we shall comment on it below. 

The most conservative approach to the Hamiltonian modification is suggested by ([3]). 
Since the problem arises due to the impossibility of integrating f(t,e) over t at e — 0, the 
natural solution is to modify f(t, e) locally around (in the e- neighborhood of) t = in such 
a way that the integral can be taken (note that we are leaving the operator structure of the 
Hamiltonian intact). 



1 This assumption is actually stronger than one might naively have thought: seemingly small interaction 
terms in the Hamiltonian are sometimes responsible for large quantum effects, for instance due to degrees 
of freedom becoming light. An example where this happens is the matrix big bang model [Sj [9] , where an 
important one-loop potential is generated in the weak coupling region of the field theory. For this reason, 
our present discussion will not directly apply to the matrix big bang model, though we hope to treat that 
model using similar techniques in future work. 
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The subtractions necessary to appropriately modify f(t,e) are familiar from the theory 
of distributions. Namely, for any function f(t, e) developing a singularity not stronger than 
l/t p as e is sent to 0, with an appropriate choice of c n (e), one can introduce a modified 

f{t,e) = f(t,e)-*jr Cn {e)5V(t) (4) 

n=0 

(where 6^(t) are derivatives of the 5-function) in such a way that the e — > limit of f(t, e) 
is defined in the sense of distributions. The latter assertion would imply that the e — ► 
limit of 

'' f(t,e)F(t)dt (5) 



is defined for any smooth "test-function" !F(t), and, in particular, that the e — > limit of 
([3]) becomes well-defined, if f(t, e) is replaced by f(t, e). (Note that, since f(t, e) and f(t, e) 
only differ in an infinitesimal neighborhood of t = 0, this modification will not affect the 
evolution at finite t). 

As a matter of fact, the subtraction needed for our particular case is simpler than (j4j). 
Since the n > terms in (J4]) can only affect the value of the evolution operator (j3J) at t' = 0, 
if one is only interested in the values of the wave function for non-zero times, one can simply 
omit the n > terms from (j3J). One can then write down the subtraction explicitly as 



/(*,£) = /(*,£)- / f(t,e)dt\6(t). (6) 




The appearance of a free numerical parameter (which can be chosen as to in the expression 
above, or a function thereof) is not surprising, since, if f(t,e) is an adequate modification 
of f(t,e), so is f(t,e) + c5(t) with any finite c. 

For the particular l/\t\ time dependence of the Hamiltonian mentioned in the introduc- 
tion, one can choose f(t,e) as l/Vt 2 + £ 2 , in which case f(t,e) becomes 

/i/| t |(t,e) = -=L= + 21n(//e)<y(f). (7) 



e 



It is sometimes more appealing to replace the <5-function in ([7]) by a resolved 5-function, in 
which case we find 

(with fj, being an arbitrary mass scale). 

One should note that it is very natural to think of the above subtraction procedure as 
renormalizing the singular time dependence of the Hamiltonian. Indeed, the mathematical 
structure behind generating distributions by means of 5-function subtractions is precisely 
the same as the one associated with subtracting local counter-terms in order to render 
conventional field theories finite. For concreteness, consider the one-loop contribution to 
the full momentum space propagator in A0 3 field theory, given by the diagram 
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If we compute it using position space Feynman rules, we find that it is proportional to the 
Fourier transform of the square of the scalar field Feynman propagator D(x,x'). However, 
while the Feynman propagator itself if a distribution, its square is not. For that rea- 
son, if one tries to evaluate the Fourier transform, one obtains infinities, since integrals of 
[D(x, x')] 2 cannot be evaluated. The problem is resolved by subtracting local counter-terms 
from the field theory Lagrangian, which, for the above diagram, would translate into adding 
S(x — x') and its derivatives (with divergent cutoff-dependent coefficients) to [D(x, x')] 2 in 
such a way as to make it a distribution. The mathematical structure of this procedure is 
precisely the same as what we employed for renormalizing the singular time dependences 
in time-dependent Hamiltonians. 

We should remark upon the general status of our Hamiltonian prescription viewed 
against the background of all possible singularity transition recipes one could devise. If 
the only restriction is that the evolution away from the singularity is given by the original 
Hamiltonian, one is left with a tremendous infinitefold ambiguity: any unitary transforma- 
tion can be inserted at t — and the predictive power is lost completely. One should look 
for additional principles in order to be able to define a meaningful notion of singularity 
transition. 

Our prescription can be viewed as a very conservative approach, since it preserves the 
operator structure of the Hamiltonian (the counter-terms added are themselves proportional 
to h, the singular part of the Hamiltonian). In the absense of further physical specification, 
this approach appears to be natural and can be viewed as a sort of "minimal subtraction" . 
However, under some circumstances, one may be willing to pursue a broader range of 
possibilities for defining the singularity transition. For example, one may demand that the 
resolution of the singular dynamics must have a geometrical interpretation (at finite values 
of e). This question will be addressed in [T4] . 

In section 3, the focus of our attention will be a particular quantum system with a 
Hamiltonian quadratic in the canonical variables. For such linear systems, it is most com- 
mon to analyze quantum dynamics in the Heisenberg picture, rather than in the Schrodinger 
picture we have employed above for the purpose of describing our general formalism. For 
convenience, we shall give a summary of the relevant derivations in appendix A. In short, 
one should construct the most general classical solution of the system in the form 

x{t) = Au{t) +A*u*{t). (9) 

The solution to the Heisenberg equations of motion is simply obtained by replacing the 
integration constants A and A* in the above expression by creation-annihilation operators 
a and a', which (with an appropriate normalization of u(t)) satisfy the standard commu- 
tation relation [a, a*] = 1. The question of solving for the quantum dynamics is then most 
commonly phrased in terms of constructing the mode functions u(t) and u*(t), which are 
normalized solutions to the classical equations of motion. 
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Our prescription may equally well be applied in such setting. One can analyze the 
classical equations of motion derived from the time-dependent Hamiltonian. It is safest to 
do so at finite e, since the naive e — > limit of the classical equations of motion may not 
necessarily exist. However, the e — ► limit of the solutions for the mode functions will 
exist, and will, of course, define the same quantum dynamics as the general solution to the 
Schrodinger equation given by (J2]). 



3 Free fields on the compactified Milne universe 

3.1 The compactified two-dimensional Milne universe 

The two-dimensional Milne universe 

ds 2 = -dt 2 + t 2 dx 2 , (10) 

with < t < +oo, corresponds to the "future" quadrant X ± > of Minkowski space 
ds 2 = —2dX + dX~ via the identification 

X± = i/2 6±X - (U) 
The Milne universe can be compactified by the identification 

x~x + 2vr, (12) 
which corresponds to the discrete boost identification 

X ± ~ e ±2n X ± . (13) 

The resulting space is a cone, which is singular at its tip t = 0. 

The action for a free scalar field in the (compactified) Milne universe is 

The corresponding equation of motion 

+ ^-^ + m 2 (p = O (15) 



is solved by 



| r i( JZLx~e~ w +—X + e w +lw , 

i e ^ e -i^ H m^ mt ) (17) 



2V2 
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and their complex conjugates [TTj [T2] . Here H^ 1 ' denotes a Hankel function, and the 
compactification ffl2|) enforces the momentum quantization condition I e Z. For solutions 
to the equation of motion ffl5|) . we define the scalar product [TT] 

r27r 

— i / <ixt 



&(t, X)0*(t, z) -&(t, X)<%(t,x)\ (18) 

and the Klein-Gordon norm (0, 0). The solutions (fT6l) are normalized to have Klein-Gordon 
norm —1. 

To quantize the scalar field 0, one expands 

<f>(t,x)= [akU k (t,x) + a\u* k (t,x) , (19) 
ke7L 

where the u k (x,t) have Klein-Gordon norm 1, which ensures the canonical commutation 
relations 

[a k ,a\] = S k ,i- (20) 

We choose 

u k (t,x) =<ip^ k (t,x). (21) 

Essentially because ipm,k of ( |T6l) are superpositions of negative frequency waves on the cover- 
ing Minkowski space, the vacuum state defined with the creation and annihilation operators 
of ([191) is an adiabatic vacuum of infinite order [TT]. Note, however, that in a compactified 
Milne universe (where globally defined inertial frames are absent) this particular adiabatic 
vacuum is no more special than any other adiabatic vacuum of infinite order (of which there 
are infinitely many). 

Near t — 0, the I ^ mode functions behave as (see, for instance, [7]) 



ui 



IJItxI sinh(7r/) 



( 'mt\ d iri • / mt\ ll m,- n 

(-) + (-) 



(22) 



with (fi defined by e lLpl = T(l + il) J 8 *y and satisfying cp_i = —(pi, while 



1 / 2i, fmtW 
«o~-^= l--log — • (23) 



2v/2 V 7r to V 2 / 

The mode functions are clearly singular at t — 0. The question we now want to address is 
whether quantum mechanical evolution can be consistently and naturally defined beyond 
t = 0. 

In the literature (see, for instance, [121 H31 this question has been addressed by 
extending the range of the t coordinate in the compactified Milne metric (fTOl) to — oo < t < 
oo, i.e. by adding a "past cone" to the "future cone"J§ In the action (fT4"|) . the factor t is 
replaced by \t\, 



2 In some string theory contexts, it is natural to consider the full Minkowski space up to the discrete 
boost identification (Tl3|) . which in addition adds "whisker" regions with closed timelike curves. We will not 
consider whisker regions in the present paper. 
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The same goes for the scalar product ffl~8]) 



—i 



2- 



dx \t\ 



X)4>l(t, X) - 01 (t, X)(j>l(t, X) 



(25) 



and the corresponding Klein-Gordon norm. 

The question then is how to define matching conditions between t < mode functions 
and t > mode functions, i.e. how to define global mode functions. Natural globally defined 
mode functions are obtained by allowing X to be either both positive or both negative in 
f TTBT) (see (T2TT) ) . As these are superpositions of negative frequency Minkowski modes, they 
describe excitations above the (adiabatic) vacuum inherited from Minkowski space. The 
solutions ([TBI) have the property that they are analytic in the lower complexified t-plane. 
For t < 0, they can be written as 



ip m ,i(t,x) 



1 



*2 -ilx {ttW 



2^2 



e 2 e 



mt 



(f <0) 



(26) 



which still has Klein-Gordon norm —1. For t approaching from below, we have for the 
corresponding mode functions Ui(t,x) = ^ z (t,x) with I 7^ 0, 



Six 



Ul 



2J2nl sinh(TrZ) 



mt 
~2~ 



+ 



mt 
~2 



and 



u 



( 



, 2z n 
1 + — log 

7T 



mt 



(27) 



(2? 



2V2 V 

Note that, even though the above prescription may seem natural, and it does define 
consistent matching conditions and a unitary evolution, it should not be given any privileged 
status. The (compactified) Milne universe contains a genuine singularity at the origin, and 
the question of how the system evolves in the neighborhood of the singularity cannot be in 
principle settled through an appeal to a flat Minkowski space (even though there is nothing 
wrong with using the covering Minkowski space for constructing particular evolutionary 
prescriptions). As we shall see below, more general rules for singularity crossing can be 
devised, with a different set of mode functions and a different vacuum state (which, being 
an adiabatic vacuum of infinite order, is no better and no worse than the one inherited from 
the covering Minkowski space). 

Even though the modefunctions u\ = ipmi constructed above solve the equations of 
motion derived from the action (1241) at all positive and all negative t, there are no meaningful 
equations of motion satisfied at t = 0. Correspondingly, even though the quantum evolution 
defined in terms of the above prescription for the mode functions is unitary (and essentially 
inherited from the covering Minkowski space), this quantum evolution cannot be represented 
as a solution to the Schrodinger equation for the Hamiltonian derived from (T2"4"|) . In what 
follows, we shall nevertheless be able to cast this quantum evolution in a Hamiltonian form 
by appropriately renormalizing the time dependences in the Hamiltonian of the system. 
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3.2 Quantum Hamiltonian evolution across the Milne singularity 



In section 2, we constructed a general prescription which allows to define a Hamiltonian 
evolution across an isolated singularity in the time dependence of the Hamiltonian. Since 
the case of a free scalar field on the Milne orbifold falls precisely into this category, it will 
be instructive to compare the above consideration in terms of the covering Minkowski space 
with our general prescription. We shall see that the two are in fact related, even though it 
is only in the parametrization of section 2 that the evolution has a manifestly Hamiltonian 
form at t = 0. 

The Hamiltonian corresponding to the action (1141) is 



H 



L>2 



+ 



m \t\ 



dx i 



(29) 



Following the general guidelines presented in section 2, we shall regulate the l/\t\ time 
dependence into fi/\t\(t,e) of (jHJ): 



H=±f m (t,e) J dx (vrJ + 0' 2 ) + ^ J dx, 



(30) 



Near the origin, where the mass term is negligible, the equations of motion take the 
form 



h/\* 
f: 



it 



AW 



o 



or, after Fourier-expanding y/27T(j)(x,t) = J24>i(t) exp(ilx) 

i fy\t\ 

0- 7 — 

h/\t\ 

The general solution to this equation is 



+ 1 /i/|t|< 



0. 



Aiexp 



il / f 1 /\ t \(t,e)dt 



+ Bi exp 



H / h/\t\(t,e)dt 



or 



A exp 



il farcsinh- 

. V e 



2 . . t 
— In(//£J arctan - 

7T £ 



B, 



, / • , t 2 , . , t 
—il arcsmn — I — m(ue) arctan - 

en e 



(31) 



(32) 



(33) 



With e explicitly taken to 0, this becomes 



Bi\2fjit\ 



-il sign(t) 



(34) 



(35) 



To construct the Heisenberg field operator (which contains all information on quantum dy- 
namics) one should choose any such complex solution and, after normalizing appropriately, 
promote it to a mode function, as in ffl9l) (see also appendix A). 

The question that will interest us here is how the quantum dynamics described by the 
Hamiltonian with our "minimal subtraction" is related to the mode function prescription 
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plj) inherited from the covering Minkowski space. To this end, we shall define mode 
functions uf that solve (1321) and coincide with ui of (I2ip for t > 0; however, they will 
generically differ from ui for i < 0. To see the relation between u\ and u/, we construct 
by choosing Ai and in (I35I) in such a way that it equals (122]) for t > and then 



(/■<) 



compare it, for t < 0, with (T2j 

In order to match (|22|) and (|35|) for t > 0, we impose 



B, 



1 



2y / 2?rZ sinh(TrZ) V 4 ^, 
1 ( m\ 



_ Tri 

e 2 ' 



2J27r/sinh(vr/) V 4 ^, 



(36) 
(37) 



Then, at t < 0, 



,,ilx 



- 






m 



e 2 



2^/27rZsinh(7r/) 

Comparing this expression with (1271) . we conclude that they are indeed the same if 



(3« 



T 



exp 



-2y; + 7T 

2/ 



(39) 



Note that the fact that \x depends on the Milne momentum / implies that the value of 
the arbitrary parameter introduced by our renormalization procedure is different for each 
of the oscillators comprising the field. For that reason, even though the covering Minkowski 
space prescription turns out to be the same as our "minimal subtraction" for each of the 
oscillators, for the entire field it is not. Phrased in the Hamiltonian language, the covering 
space prescription for the Milne singularity transition turns out to be different from the 
simplest consistent recipe one could devise, even though it is related to such simple recipe 
in a fairly straightforward way. 



4 Conclusions 



We have addressed the issue of how one can define a unitary quantum evolution in the 
presence of isolated singularities in the time dependence of a quantum Hamiltonian. If 
one demands that the operator structure of the Hamiltonian should be unaffected by regu- 
larization prescriptions (the "minimal subtraction" recipe), one discovers a one-parameter 
family of distinct quantum evolutions across the singularity. 

For the case of free quantum fields on the Milne orbifold, the covering Minkowski space 
considerations previously brought up in the literature [TH [TJl [131 El t urn ou t t° be closely 
related to, though distinct from, our "minimal subtraction" proposal. One explicit advan- 
tage of our present approach is that it makes the evolution across the singularity manifestly 
Hamiltonian, which was not the case in the context of the previous discussions. 
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A Linear quantum systems 

In this appendix, we shall review the dynamics of linear quantum systems. This material is 
very basic and well-known; however, it is usually presented in relation to a few specific linear 
systems of physical interest, whereas, for our purposes, it shall be convenient to summarize 
here the treatment of a general one- dimensional linear quantum system described by the 
Hamiltonian 

H = Mp 2 + ^-X 2 (40) 

with f(t) and g(t) being arbitrary functions of time. 
The equations of motion take the form 

P = -g(t)X X = f(t)P (41) 

or 

■V • j .V • ///.V 0. (42) 

Should one succeed finding a complex solution u{t) to this equation, one would be able to 
write down the most general real solution in the form 

X(t) = Au(t) + A*u*(t) (43) 

with some complex constant A. In the quantum case, the solution to the Heisenberg 
equations of motion will have the exact same form with A and A* replaced by Hermitean- 
conjugate operators a and at: 

X H (t) = au(t) + a f u*{t). (44) 

Our solution for the quantum dynamics shall be complete if we establish the commuta- 
tion relations for a and aJ. Before doing so, we recall the important notion of Wronskian 
for a linear differential equation. For any two solutions xi(t) and x^it) of a second order 
differential equation, their Wronskian is defined as 



W[x l (t),x 2 (t)} = det 



Xi x 2 

Xi X 2 



(45) 



It is straightforward to show that, for equation ( 1421) . the Wronskian of any two given 
solutions satisfies 

W = jW. (46) 
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In other words, W/f does not depend on time. This circumstance permits to define the 
"Wronskian norm" for any complex solution u(t): 



IML = -iBf3. ,47) 

As we have just demonstrated, the value of this expression does not depend on the moment 
of time one chooses to evaluate it. The familiar Klein-Gordon norm for free quantum fields, 
which we use in section 3, is a direct generalization of the Wronskian norm. 

The physical relevance of the Wronskian norm becomes apparent from the consideration 
of commutators: 

j Till* uii* 

-i = [P(t),X(t)} = j{ X ,X\ = J K flt ] = -iNL[a,a + ]- ( 4 §) 

Therefore, to obtain the standard commutation relations for the creation-annihilation op- 
erators, [a, a'] = 1, one has to choose a complex solution u(t) with Wronskian norm 1. 
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